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ABSTRACT 


Eight  earthquakes  from  the  Caucasus  occurring  from  1971  to  1975  were  ex- 
amined in  a seismic  discrimination  context.  Seismograms  from  ALPA,  LASA, 

NORSAR,  and  the  HCLP  and  the  WWSSN  stations  were  studied  for  source  mechanism, 
corner  frequency,  pP,  complexity,  and  spectral  ratio.  All  the  Caucasus 
events  can  be  identified  as  earthquakes  by  means  of  pP  and  Mg:m^  when  compared  to 
Kazakh  and  southwest  USSR  explosions.  However  it  seems  plausible  that  these  dis- 
criminants could  be  rendered  ineffective  for  events  in  the  Caucasus  by  means  of 
a shot  array. 
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INTRODUCTION 


Discrimination  parameters  have  been  applied  in  detail  to  events  from 
only  a few  nuclear  test  sites  and  earthquake  regions.  This  study  is  part  of 
a series  which  extends  these  discrimination  studies  to  other  regions  of  shal- 
low earthquakes.  This  particular  report  examines  eight  such  earthquakes  in 
the  Caucasus  of  the  USSR.  Most  of  the  earthquakes  in  the  Caucasus  are  prob- 
ably associated  with  large  thrust  faults,  which  are  a result  of  the  Meso- 
Cenozoic  Alpine  orogeny.  Actual  seismogram  analysis  revealed  that  the  events 
chosen  were  too  small  to  determine  their  fault  plane  solutions.  We  determine 

average  M and  m,  for  the  events  in  this  studv  and  apply  other  common  dis- 
s b 

crimination  parameters  such  as  first  motions,  corner  frequency,  complexity, 
and  short-period  spectral  ratios  ft)  the  data.  Results  for  the  Caucasus  events 
will  be  compared  to  Kazakh  and  USSR  explosions. 


TECTONIC  SETT INC 


Coiu-ral  Features 

Figure  1 is  a map  of  the  Caucasus  region.  The  principal  geologic  features 
in  this  area  trend  northwest-southeast  and  are  associated  with  the  Meso-Cenozoic 
Alpine  orogeny.  In  general  the  seismicity  follows  the  trend  of  the  mountains, 
hut  zones  of  higher  activity  occur  near  AIN,  A7K  and  AIN,  A8E  (Nowroozi,  1*171). 

In  general  the  seismicity  is  well  enough  distributed  that  no  area  is  exempt  from 
high  seismic  risk  according  to  the  Russian  observations  summarized  in  Byus  et  al. 
(19b8).  The  majority  of  the  earthquakes  occur  within  the  crust  in  the  Caucasus 
Mountains  (Keil is-Borok,  I960;  Byus  et  al.,  1968;  Nowroozi,  1971);  but  where  the 
seismic  zone  extends  into  the  Caspian  Sea,  the  earthquakes  are  subcrustal,  with 
depths  ranging  from  50  to  150  km.  These  deeper  events  are  possibly  due  to  the 
subsidence  of  the  Caspian  Basin  which  has  accumulated  as  much  as  20  km  of  sedi- 
ment in  its  southern  portion. 

Source  Mechanisms  of  Earthquakes 

NE1S  epicenters  for  events  in  this  area  from  January  1971  through  February 
1975  are  shown  in  Figure  1.  The  geologic  faults  in  this  figure  were  inferred 
from  a satellite  photomosaic  of  the  area.  Most  earthquakes  are  associated  with 
thrust  faults  in  this  fold  belt  since  compress iona 1 axes  determined  from  fault- 
plane  solution  are  generally  perpendicular  to  the  trend  of  the  geologic  struc- 
ture (Keil is-Borok,  1960;  Nowroozi,  1971).  Within  this  dominant  NW-SE  strike, 
there  is  considerable  variance  to  the  fault-plane  orientations;  and  strike-slip 
motion  is  not  uncommon. 

Velocity  Model 

Deep  seismic  sounding  studies  report  crustal  thicknesses  reaching  50  to  60 


Nowroozi,  A.,  1971.  Seismo— tectonics  of  the  Persian  Plateau,  Eastern  Turkey 
Caucasus,  and  Hindu-Kush  Regions,  Bu 1 1 . Se i sm.  Soc.  Am. , 61,  317. 

Byus,  E.  1.,  A.  D.  Tskhakaya,  and  M.  M.  Rubinshtein,  19b8.  Regional  seismic 
zoning  of  the  USSR:  Georgia,  in  Seismic  Zoning  of  the  USSR,  edited  hv 
S.  V.  Medvedev,  English  translation,  Kotor  Publishing  House,  Jerusalem, 
Israel. 

Koi 1 is-Borok,  V.  L.,  I960.  Investigation  of  the  Mechanism  of  Earthquakes, 
English  translation,  American  Geophysical  Union,  Washington,  D.  G. 


USSR 


km  in  the  Caucasus  (BeUayevsky  et  al.,  1968;  Borusk  and  Shalpo,  1976)  and  35 
to  40  km  in  the  central  Caspian  Sea  (Neprochnov,  1968).  Table  I shows  the  gen- 
eral crustal  and  upper  mantle  velocity  structure  for  the  Caucasus  area  of  Fig- 
ure 1 (Beliayevsky  et  al.,  1968). 
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Earthquake  Selection 

For  this  study  shallow-focus  earthquakes  were  selecteJ  from  the  Caucasus 
Mountains  and  Caspian  Sea  regions.  The  NEIS  list  suggests  that  the  large  maj- 
ority of  the  earthquakes  in  the  Caucasus  mountains  occur  within  the  crust. 

Where  the  mountains  extend  into  the  Caspian  Sea,  .he  earthquakes  are  subcrus- 
tal,  with  depths  ranging  from  50  to  150  km.  The  earthquakes  chosen  were  lim- 
ited to  the  years  1971  to  1975  so  that  we  could  utilize  the  data  from  the  large 
seismic  arrays  and  the  HGLP  network.  Initially,  earthquakes  with  low  M^, 
relative  to  m^,  were  sought;  but  a survey  of  available  data  revealed  no 
earthquakes  with  anomalously  low  in  this  area  for  1971-1975.  Those  earth- 
quakes with  the  lowest  observed  values,  again  relative  to  m^,  (NORSAR 
observations)  were  simply  selected  out.  We  selected  a total  of  8 earthquakes 
as  listed  in  Table  II  and  shown  in  Figure  I.  Table  II  also  lists  the  depths 
determined  from  pP  observations  at  LASA,  NORSAR,  and  selected  WWSSN  stations. 

All  of  the  selected  events  are  located  within  the  crust,  except  the  Caspian 
Sea  earthquake  (event  5). 

Explosion  Selection 

The  Caucasus  earthquakes  were  compared  to  explosions  in  Kazakh  and  south- 
west USSR  which  lie  to  the  north  and  northeast  of  the  Caucasus.  We  selected  a 
total  of  8 explosions  as  listed  in  Table  III  for  this  comparison. 

Seismic  Stations 

Digital  data  from  the  three  arrays  ALPA,  LASA,  and  NORSAR  and  from  the 
available  HGLP  stations  and  film  data  from  selected  WWSSN  stations  were  assem- 
bled for  these  events.  WWSSN  stations  were  selected  on  the  basis  of  magnifica- 
tion and  proximity  to  the  region  of  study,  and  it  is  unlikely  that  stations  not 
used  here  would  significantly  add  to  the  data  base  for  these  Caucasus  events. 

The  available  I.ASA  and  NORSAR  short-period  P recordings  are  shown  in  Figure 
2.  We  have  used  the  AO  subarray  from  LASA  and  the  C3  subarray  from  NORSAR  as  in 
our  previous  reports  of  this  series  in  order  to  avoid  signal  loss  at  high  fre- 
quencies due  to  fuil-array  beaming.  The  subarray  beams  have  individual  sen- 
sors delayed  according  to  the  appropriate  velocity. 
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TABLE  II 


NEIS  Parameters  for  Caucasus  Earthquakes 


Origin  Time  LAT  (N8) 


7106280 

710908 

711015 

720203 

731214 

741223 

750109 

750220 


10:53:45.7 

22:35:15.8 

17:08:06.3 

02:29:21.9 

09:11:46.3 

05:22:08.4 

23:09:46.6 

14:44:25.6 


42.400 

41.100 

41.428 

40.719 

41.876 

43.128 

42.889 

42.431 


lonc;  ( ° E ) 

43.300 

43.800 

48.571 

48.403 

49.028 

47.049 

46.987 

45.179 


Depth 

This 

NEIS  Stuch 


TABLE  III 


NEIS  Parameters  for  Kazakh  and  southwest  USSR  Explosions 


Event 


Date 

Origin  Time 

LAT  ( °N) 

LONC  ( °E) 

m. 

711222 

06 

59 

56.3 

47.872 

48.222 

— n 

6.0 

721003 

08 

59 

57.8 

46.848 

45.010 

5.8 

721124 

09 

59 

57.8 

51.843 

64.152 

5.2 

730815 

01 

59 

57.8 

42. 711 

67.410 

5.3 

730919 

02 

59 

57.2 

45.635 

67.850 

5.2 

730930 

04 

59 

57.5 

51.608 

54.582 

5.2 

740625 

03 

56 

57.6 

49.889 

78.115 

4.8 

741207 

05 

59 

56.9 

49.908 

77.648 

4.8 

SICNAI.  ANALYSIS 


Source  F.ttocts 

We  wanted  initially  to  Identify  a source  mechanism  for  each  event  so  that 
we  could  predict  the  radiation  pattern  for  body-wave  and  surface-wave  phases. 
Determination  of  the  fault  planes  for  earthquakes  whose  magnitudes  are  below 
b.l)  Is  generally  not  reliable  with  teleseismlc  data,  and  all  our  earthquakes 
have  nij}  less  than  b.  first  motions  were  determined  when  possible  from  the 
data  we  had  collected,  and  International  Sc i sinolog  leal  Center  (ISC)  bulletin 
data  was  available  for  events  l through  5.  Wo  have  plotted  all  the  first 
motions  lor  the  Caucasus  earthquakes  in  figure  3.  Unfortunately,  the  short- 
period  data  were  not  consistent,  and  it  is  impossible  on  the  basis  of  the 
plots  alone  to  determine  iiul  i v (dun  1 t an 1 1 planes  because  dilatational  and 
compress  Iona  1 tirst  motions  do  not  separate  well.  However,  taken  together, 
there  Is  a dominance  ol  compress  Iona l motion  in  the  central  to  NK  central 
part  ol  the  local  sphere  plots;  this  would  result  from  thrust  faulting  NK  in 
a direction  perpendicular  to  the  trend  ot  the  Caucasus  Mountains  and  is  con- 
sonant with  the  body  of  known  focal  mechanisms  tot-  this  region. 

Corner  frequencies  and  seismic  moments  for  our  Caucasus  events  have  been 

estimated  from  short-period  1J\SA  and  NORSAR  spectra.  The  spectra  are  shown  in 

Figure  4.  These  spectra  are  from  the  phased  beams  of  the  AO  subarray  at  l.ASA 

and  Cl  subarray  at  NOKSAK,  The  sample  length  was  t>.4  seconds.  The  signals 

have  been  tapered,  and  the  Instrument  response  was  removed  from  the  signal  and 

noise  spectra;  but  noise  spectra  have  not  been  subtracted  from  the  signal 

spectra.  Attenuation  was  removed  from  the  observed  (t*  «*  0.0)  spectra  by 

_ ) 

multiplying  by  the  factor  exploit*)  with  a t*  ol  .b)  for  a io  *"  source  model 

-1  -2  -3 

and  .44  lor  a u>  source  model  at  l.ASA  and  t*  ot  .28  lor  w and  .09  lor  w 

source  models  at  NOKSAK.  The  basis  lor  these  t*  values  will  be  shown  Later 

in  this  report.  Corner  frequencies  were  estimated  with  the  assumption  ol 

-2  -3 

complete  stress  drop  and  a w or  w asymptotic  relation  at  high  frequencies. 
Many  of  the  estimates  here  ot  corner  frequency  and  long-period  level  are  ten- 
uous, and  those  which  were  made  with  little  confidence  are  represented  by 
dashed  lines  in  Figure  4.  Occasionally  we  computed  spectra  lor  25b  second 
windows  for  those  cases  where  the  si gna 1-to-noise  ratio  was  adequate  at  lower 
frequencies.  For  these  cases  the  shape  ot  the  spectra  was  sufficiently  simi- 
lar to  those  for  the  b.4  second  windows  that  wo  have  not  shown  them. 
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event  1 


Hrst  motions  for  Caucasus  earthquakes  in  the  I 
focal  sphere  (Wulff  net). 


over  halt  of  the 


EVENT  2 


Figure  ) (cont.) 


1'lrst  motions  for  Caucasus  earthquakes  in  the  lower  half  of 
the  focal  sphere  (Wulff  net). 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
. DISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


FREQUENCY  (Hz) 

Figure  4.  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Caucasus 
earthquakes  with  instrument  response  and  attenuation  removed. 


FREQUENCY  (Hi) 

SA  AO  and  NOKSAR  C3  subarray  spectra  of  P waves  from  Caucasus 
rthquakes  with  instrument  response  and  attenuation  removed. 
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LASA  AO  and  NORSAK  C3  subarray  spectra  of  P waves  from  Caucasus 
earthquakes  with  instrument  response  and  attenuation  removed. 
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Figure  4 (cont.)  LASA  AO  and  N'ORSAR  C 3 subarray  spectra  of  P waves  from  Caucasus 
earthquakes  with  instrument  response  and  attenuation  removed. 
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SA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Caucasus 
rthquakes  with  instrument  response  and  attenuation  removed. 
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Figure  4 (cont.)  LASA  AO  and  NORSAK  C3  subarray  spectra  of  P waves  from  Caucasus 
earthquakes  with  instrument  response  and  attenuation  removed. 
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Figure  4 (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Caucasus 
earthquakes  with  instrument  response  and  attenuation  removed. 


Seismic  moment  was  calculated  from  the  long-period  body-wave  spectral 

displacement  !u  I estimated  in  Figure  4 using  the  relation 
o 
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where  R4  was  assumed  to  be  unity  due  to  lack  of  knowledge  of  the  focal  mech- 
anism, values  of  pi’  from  Table  1 appropriate  to  the  source  depth  were  used, 
and  the  divergence  factor  D was  applied  (Ben-Menahem  et  al.,  1965).  A graph 
of  moment  versus  corner  frequency  is  shown  in  Figure  5.  Our  events  fall  close 
to  or  between  the  constant  stress  drop  lines  of  Hanks  and  Thatcher  (1972), 
indicating  that  our  events  are  of  intermediate  stress  drop. 

We  attempted  to  find  the  LK  radiation  pattern  for  each  event  by  plotting 
the  anti  logs  of  the  M values  for  each  event  as  listed  in  Table  IV.  The 
results  for  the  largest  M earthquake,  event  7,  are  shown  in  Figure  6.  The 
maximum  antilog  is  plotted  on  the  outer  circle,  and  all  the  other  antilog 
values  are  scaled  linearly  between  the  center  of  the  plot  and  the  radius  of 
the  circle.  The  observed  distribution  of  amplitudes  does  not  obviously  con- 
form to  any  quadripole  radiation  pattern,  which  is  probably  due  to  varying 
propagation  effects  coupled  with  the  large  epicentral  distances  and  the  small 
magnitude  of  the  events  studied.  The  scatter  of  data  points  for  all  the  events 
was  so  poor  that  no  attempt  was  made  to  find  the  best-fitting  radiation  pat- 
terns. 

Propagation  Kffects 

We  measured  relative  mantle  attenuation  along  the  paths  from  the  Caucasus 
to  LASA  and  NORSAR,  the  only  two  sites  for  which  we  had  short-period  digital 
data.  Our  measurement  of  attenuation  is  the  factor  t*,  which  is  the  travel 
time  along  the  ray  path  divided  by  the  average  attenuation  factor  Q. 


Ben-Menahem,  A.,  S.  Smith,  and  T.  Teng,  1965.  A procedure  for  source  studies 

from  spectrums  of  long-period  seismic  body  waves.  Bull.  Seism.  Soc . Amer . , 
55,  203. 

Hanks,  T. , and  W.  Thatcher , 1972.  A graphical  representation  of  seismic  source 
parameters,  .J.  C.cophvs , Res. , 77  , 4393. 
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TABLE  IV 


Magnitude  Data  for  Caucasus  Earthquakes 


Event  I 


Stat Ion 

A 

mb 

M 

s 

HD1. 

64.0 

4.76 

cue 

52.3 

4.21 

COL 

72.3 

4.44 

USA 

87.2 

4.56 

MAT 

69.4 

5 . 50 

NAI 

44.1 

4.  19 

Nil. 

25.0 

4.56 

NORSAH 

26.9 

4.19 

NUR 

21.3 

4.10 

que 

22.5 

3.78 

SDH 

63.5 

4.68 

Sill. 

43.0 

4.16 

TOL 

35.3 

4.85 

4.44 

AV CRACK 

Event  2 

4.46 

Stat ion 

A 

mb 

M 

s 

HAG 

69.5 

4.91 

BUI. 

62.8 

4.97 

cm; 

51.8 

4.76 

COL 

73.6 

4.37 

EIL 

13.6 

4.29 

ROD 

42.8 

4.63 

RON 

28.0 

4.57 

4.17 

USA 

88.6 

3.93 

MAI. 

37.3 

4.17 

NIL 

24.4 

5.33 

NOKSAR 

28.2 

4.88 

3.62 

NUR 

22.7 

4.53 

OGD 

80.  7 

4.44 

POO 

34.2 

4.76 

3.94 

QUE 

21.7 

4.05 

3.  79 

SHI. 

42.5 

4.56 

TOL 

35.9 

4.70 

3.84 

AVERAGE 

4.65 

1.92 
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Event  3 


Stat ion 

A 

AI.I’A 

7 1.0 

II  AC 

65.8 

IUII. 

64.  1 

CMC 

48.3 

COI. 

72.  <) 

Ell. 

16.2 

FBK 

73.0 

CDH 

Vi.  7 

KOI) 

40.2 

KON 

30.0 

IJ\SA 

89.4 

MAI. 

40.8 

MAT 

66.5 

Nil. 

21.0 

NORSAR 

30.0 

NIIR 

24.0 

(XU) 

82.8 

POO 

31.4 

QUE 

18.6 

SHI. 

38.9 

TLO 

AVERAGE 

39.3 

5.00 

4 . 66 

3.67 

4 . 99 

5.04 

3.66 

3.80 

5.25 

4.42 

2.96 

3.88 

4.  /z 

3.  73 

5.28 

4.84 

3.51 

4.84 

5.12 

4.83 

4.75 

4.87 

4.58 

3.58 

4.45 

3.93 

3.42 

4.91 

3.63 

3.09 

4.76 

3.71 

■ 1 


I 

> 


Kvent  4 


Stilt  Ion 

A 

AI.PA 

73.7 

Al.Q 

100.9 

Al.Q 

100.9 

RAG 

66.0 

BUI. 

63.6 

CHC 

48.2 

CHG 

48.2 

COI. 

73.6 

(7TA 

108.7 

Ell. 

15.6 

KBK 

74.0 

KON 

30  5 

KON 

30.5 

LASA 

90.  1 

MAT 

67.0 

NAI 

43.3 

Nil. 

20.9 

NOR 

46.2 

NORSAR 

30.5 

NHR 

24.6 

OGI) 

83.3 

QUE 

18.3 

SOB 

64.3 

SHI. 

38.9 

TOI. 

39.4 

AVERAGE 

4.67 

5.22 

5.43 

5.00 

4.54 

4.96 

4 . 46 

4.60 

5.31 

4.94 

5.22 

4.47 

4.60 

4.11 

5.09 

5.67 

4.57 

5.38 

4.68 

4.  70 

4.79 

4.58 

5.29 

5.22 

5. 00 

5.03 

4.28 

4.45 

4.84 

4.99 

4.47 

4.88 

4.29 

4.93 

4.77 

m 
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Event  5 


Station 

A 

% 

M 

s 

ALPA 

72.5 

3.83 

ALQ 

100.0 

4.46 

BAG 

65.5 

5.26 

BUL 

64.8 

4.80 

CHG 

48.0 

4.80 

COL 

72.4 

5.16 

CTA 

108.5 

4.00 

DAG 

44.4 

5.09 

KOD 

40.3 

4.10 

KON 

29.9 

3.83 

LASA 

89.1 

5.56 

3.90 

MAT 

66.0 

4.25 

NAI 

44.5 

4.00 

NIL 

20.7 

4.94 

NORSAR 

29.9 

5.49 

3.38 

NUR 

23.8 

5.27 

QUE 

18.5 

4.01 

SDB 

65.5 

4.96 

SHL 

38.7 

4.51 

TLO 

39.6 

3.72 

AVERAGE 

4.89 

3.86 

Event  6 


Station 

A 

"b 

M 

s 

ALQ 

98.3 

4.60 

4.77 

BAG 

66.9 

5.00 

4.55 

BUL 

65.5 

5.06 

4.61 

CHG 

49.7 

4.97 

4.24 

CHG 

49.7 

4.48 

CTA 

110.1 

4.53 

DAG 

42.8 

4.60 

DAV 

76.8 

4.87 

KIP 

111.4 

4.71 

MAL 

39.4 

4.80 

4.59 

MAT 

66.7 

4.96 

4.97 

MAT 

66.7 

5.10 

NAI 

45.3 

4.87 

4.44 

NIL 

22.5 

4.79 

NUR 

22.0 

5.23 

OGD 

80.8 

4.70 

QUE 

20.5 

3.79 

SDB 

65.6 

4.57 

TLO 

37.8 

4.14 

AVERAGE 

4.86 

4.57 
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Event  7 


Station 

A 

% 

M 

s 

ALPA 

71.6 

4.78 

ALQ 

98.5 

5.07 

4.88 

ALQ 

98.5 

4.82 

BAG 

67.0 

5.54 

4.85 

BUL 

65.3 

5.00 

4.83 

CHG 

49.7 

5.30 

COL 

71.6 

5.03 

4.99 

CTA 

110.1 

4.61 

DAG 

43.0 

4.70 

GDH 

53.9 

5.09 

LASA 

87.7 

5.89 

MAL 

39.4 

4.76 

4.28 

MAT 

66.9 

5.57 

MAT 

66.9 

5.42 

NAI 

45.1 

5.01 

4.76 

NIL 

22.5 

5.29 

4.61 

NORSAR 

28.1 

5.01 

4.70 

NUR 

22.2 

5.22 

QUE 

20.4 

4.01 

SDB 

65.4 

4.96 

5.31 

TLO 

37.8 

4.14 

AVERAGE 

5.13 

4.78 

Event  8 


Station 

A 

“b 

M 

s 

ALPA 

72.2 

3.69 

ALQ 

98.4 

4.60 

4.62 

BUL 

64.4 

4.36 

CHG 

50.9 

4.54 

3.39 

COL 

72.2 

4.73 

DAG 

43.0 

4.80 

GDH 

53.8 

4.91 

KIP 

112.6 

4.58 

LASA 

87.7 

5.19 

4.01 

MAT 

68.2 

4.25 

NIL 

23.6 

4.50 

NORSAR 

27.7 

4.49 

NUR 

21.9 

4.79 

OGD 

80.4 

4.74 

4.07 

QUE 

21.3 

4.67 

3.59 

SDB 

64.3 

4.70 

SHL 

TLO 

41.6 

36.7 

4.33 

3.68 

T 7R 

AVERAGE 

4.67 

3.97 

-1*2- 


We  determine  t*  by  first  normalizing  and  then  summing  the  spectra  shown 
in  Figure  4 for  each  of  the  two  paths  - Caucasus  to  LASA  and  Caucasus  to  NOR- 
SAR.  These  averages  are  shown  in  Figure  7 and  8 and  represent  the  average 

relative  ground  displacement  density  spectra  at  the  true  sites.  If  we  assume 

_2 

that  the  source  spectra  has  an  f fall  off  at  high  frequencies,  then 

~ c e~ 2 “'"ft* 

F(f ) - S*  f • e 

for  frequenceis  above  the  corner  frequency  where  F(f)  is  the  average  of  the 
observed  spectrum  (corrected  for  seismograph  response),  S is  the  source  spec- 
trum scale  constant,  and  f is  frequency  in  Hz.  Then  after  taking  the  natural 
log  of  this  equation 

In  [F(f ) ] + 2-in(f)  = -irt*f  + in(S) 

If  we  plot  Jin  [F(f)]  + 2«£n(f)  versus  frequency  for  each  path,  as  shown  in  Fig- 
ures 9 and  11,  the  slope  of  the  graph  is  -nt*.  We  calculated  the  slope  in  the 
frequency  range  1.0  to  2.5  Hz,  where  the  signal  to  noise  ratio  was  highest, 

and  presumably  beyond  the  corner  frequency  so  that  the  high-frequency  asyptotic 
-2 

slope  of  f characterizes  the  source  spectrum.  Figure  5 shows  that  some  of 
the  corner  frequencies  are  above  1 Hz,  so  the  data  may  not  be  well  fitted  by 
the  above  equation;  but  if  a higher  frequency  range  were  used,  we  would  ser- 
iously contaminate  the  average  spectra  with  background  noise.  The  effect  of 
using  spectral  estimates  at  frequencies  below  the  corner  frequency  is  to  bias 
the  slope  or  t*  toward  a lower  value  than  the  actual  one.  Thus  our  t*  results 
here  are  not  as  reliable  as  needed  to  confidently  quantify  the  attenuation 

-3 

along  the  two  paths.  Similarly,  the  equation  for  a source  spectra  with  an  f 
fall  off  is 


in  [F (f ) ] + 3*in(f ) = -nt*f  + in(S) 

If  we  plot  the  left  hand  side  of  this  relation  versus  frequency  for  each  path 

as  shown  in  Figures  10  and  12,  the  slope  of  the  graph  is  again  -irt*;  but  it  is 

-3  -2 

based  on  an  assumed  f source  spectrum.  The  t*  values  for  an  f source  model 

are  .63  + .08  for  the  Caucasus  to  LASA  path  and  .28  + .06  for  the  Caucasus 

-3 

to  NORSAR  path.  The  t*  values  for  an  f source  model  are  .44  + .08  for  the 

Caucasus  to  LASA  path  and  .09  + .05  for  the  Caucasus  to  NORSAR  path.  We  add 

-2 

that  there  is  a predominance  of  theoretical  models  which  imply  a f source 

_3 

spectrum  rather  than  a f one  (Geller,  1976). 


Ce Her,  R.  J.,  1976.  Scaling  relations  for  earthquake  source  parameters  and 
magnitude.  Bull.  Seism.  Soc.  Am. , 66,  1501. 
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There  appears  to  be  is  a real  difference  in  attenuation  between  the  two 
paths  as  revealed  in  the  spectra;  note  that  the  waveforms  in  Figure  2 do  not 
clearly  reflect  this.  rhe  Caucasus  to  NORSAR  path  has  less  attenuation,  as 

expected  since  the  path  is  mostly  through  continental  shield  areas.  Values 

^ •) 

of  t*  for  the  explosions  for  an  f “ source  model  were  also  calculated  (von 
Seggern  and  Blandfotd,  1972).  These  t*  values  are  .35  + .17  for  the  southwest 
USSR  (events  1 and  2)  to  I.ASA  path,  .29  + .13  for  the  West  Kazakh  (events  3 
through  6)  to  1.ASA  path,  .19  + .04  for  the  West  Kazakh  to  NORSAR  path,  .16  + 
.10  for  the  East  Kazakh  (events  V and  8)  to  LASA  path  and  .05  + .08  for  the 
East  Kazakh  to  NORSAR  path.  NORSAR  signals  were  not  recovered  for  the  south- 
west USSR  explosions.  As  expected,  the  explosion  t*  values  for  the  paths  to 
NORSAR  are  lower  than  the  values  for  the  paths  to  LASA  since  the  path  to 
NORSAR  is  through  continental  shield  areas  while  that  to  LASA  passes  through 
the  more  complicated  tectonic  province  of  the  western  United  States,  where 
attenuation  is  known  to  be  high. 


von  Seggern  and  Sobel  (1975)  have  presented  a method  of  predicting  20- 
second  Rayleigh-wave  amplitudes  by  geometrical  ray  tracing.  Their  work  shows 
that  lateral  changes  in  20-second  LR-phase  velocity  can  cause  major  multipath- 
ing  and  also  large  teleseismic  amplitude  variations  due  to  accumulated  small- 
scale  refraction  effects.  Figure  13  illustrates  the  paths  of  rays  which  result 
for  event  7 (42. 9N,  47. OF).  (The  rays  which  pass  close  to  the  North  Pole  were 
stopped  close  to  that  point  due  to  ray-tracing  numerical  problems  near  the  poles.) 
No  attempt  is  made  to  quantitivelv  relate  these  patterns  to  observed  LR  amplitudes 
since  we  cannot  correct  tire  predicted  amplitude  for  source  mechanism,  but  we  can 
qualitatively  compare  the  pattern  of  rays  to  M residuals  calculated  from  the 
magnitude  data  in  Table  IV.  Where  the  rays  cross  we  can  expect  either  construc- 
tive interference  (large  positive  residuals  as  at  ALQ , MAT,  001)  and  SOB  in 
Table  IV)  or  destructive  interference  (large  negative  residuals  as  at  NORSAR, 

QUE  and  T0L) . The  M residuals  at  CMC  and  LASA  are  negative,  as  can  be  pre- 
dicted from  the  diverging  raypaths  near  those  sites.  Stations  which  record  a 
stable  amplitude  are  desirable  for  calibrating  a source  region  in  magnitude 


von  Seggern,  P.,  and  R.  Blandford,  1972.  Source  time  functions  and  spectra 
for  underground  nuclear  explosions,  Oeophvs,  J.. , 31,  83-97. 


von  Seggern,  I).,  and  P.  Sobel,  1975.  Experiments  in  refining  Ms  estimates  for 
seismic  events,  SDAC  Report  No:  TR-75-17,  Teledyne  Geotech,  Alexandria, 
Virginia. 
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determination.  I'o r events  from  tlie  Caucasus,  such  stations  could  be  located 
in  eastern  China  and  USSR  and  central  Kurope.  The  complicated  LR  ray-tracing 
t igure  partially  explains  why  we  were  unable  to  determine  LR  radiation  patterns 
from  our  amplitude  observations. 
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DISCRIMINATION  ASPECTS 


V\ 


In  this  and  the  following  sections  we  apply  several  common  discriminants 
between  earthquakes  and  explosions  to  the  Caucasus  earthquakes  and  Kazakh  and 
Southwest  Russia  explosions.  We  independently  determined  average  Mg  and  m^ 
values  for  the  Caucasus  events  in  Table  II  using  WWSSN  stations,  HOLP  stations, 
and  the  ALPA,  LASA,  and  NORSAR  arrays.  The  data  used  in  the  average  estimates 
is  listed  in  Table  IV.  Negative  values  in  the  table  indicate  noise  measure- 
ments which  were  used  as  upper  bounds  for  the  signal  amplitudes.  Magnitudes 
were  computed  according  to  the  formulas 
= log  (A/T)  + B(A) 

log  (A/T)  + 1.66  log  A + 0.3 


M 


where  A one-half  the  peak-to-peak  maximum  recorded  amplitude  reduced  to  mvi 
ground  displacement. 


T = period  in  seconds  (restricted  to  17  to  23  sec  for  Mg  calculations), 


A = epicentral  distance  in  degrees, 

B(A)  = Gutenberg-Richter  correction  term  for  surface-focus  P waves. 


We  also  independently  determined  Mg  for  the  Kazakh  and  Southwest  USSR  explo- 


sions in  Table  III  using  the  NORSAR  array  and  HGLP  stations.  The  data  used 


in  the  average  estimates  for  explosion  Mg  are  listed  in  Table  V,  and  explo- 


sion values  are  from  the  NEIS  list.  Throughout  we  have  used  a method  of 
magnitude  averaging  proposed  by  Ringdal  (1976)  in  which  the  magnitudes  at  the 
individual  stations  are  assumed  to  follow  a Gaussian  distribution.  Among  this 
distribution,  amplitudes  corresponding  to  some  magnitudes  will  fall  below  the 
noise  level;  and  Ringdal's  method  then  substitutes  a noise  measurement  at  those 
stations  which  don't  detect  and  computes  the  maximum  likelihood  estimate  of 
magnitude  based  on  measured  signals  and  noise.  The  effect  of  this  procedure  is 
to  more  accurately  define  the  magnitude  of  events  not  widely  recorded;  more 
specifically,  it  usually  results  in  a lower  magnitude  than  that  produced  by 
routine  averaging  of  observed  magnitudes.  The  Ms-m^  plot  using  this  method  for 
all  the  events  of  this  study  is  shown  in  Figure  14.  Note  that  two  explosions 
(events  7 and  8 in  East  Kazakh)  had  no  LR  detections,  and  an  upper  bound  has 


been  plotted  for  Mg  using  the  NORSAR  noise  level  of  Table  V. 


Ringdal,  F.,  1976,  Maximum- likelihood  estimate  of  event  magnitude,  Bull. 
Seism.  Soc.  Am.,  66,  789. 


-54- 


TABLE  V 


Magnitude  Data  for  Kazakh  and  southwest  USSR  Explosions 

Event  1 


Station 

A 

M 

s 

CHG 

50.1 

4.25 

EIL 

21.0 

4.25 

KON 

25.2 

4.25 

NAO 

35.9 

4.92 

AVERAGE 

4.44 

Event  2 


Station  A M 

s 

KON  24.3  3.25 

NAO  25.4  3.55 

AVERAGE  3.40 


Event  3 


Station  A M 

s 

CHG  42.8  -3.44 

KON  30.7  2.86 

AVERAGE  2.86 


Event  4 


Station  A M 

s 

CHG  35.8  -3.83 

KON  38.4  3.67 

MAT  53.3  -3.85 

AVERAGE  1 An 


Event  5 


Station 

A 

M 

s 

KON 

36.6 

3.31 

MAT 

52.1 

-4.08 

AVERAGE 

Event  6 

3.31 

Station 

A 

M 

s 

KON 

26.0 

3.54 

AVERAGE 

Event  7 

3.54 

Station 

A 

M 

s 

NAO 

38.1 

-2.61 

AVERAGE 

Event  8 
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The  magnitude  averages  of  those  events  where  many  of  the  readings  are 
noise  levels  are  lower  than  what  would  result  if  only  measured  signal  ampli- 
tudes were  used,  as  expected  from  Ringdal's  method.  Where  only  noise  measure- 
ments were  used  in  the  determination,  an  arrow  in  the  figure  represents 

the  maximum  M value.  The  line  M = m -1.5  (suggested  as  a discriminant  line 
s s b 

by  Blandford  and  Clark  (1975)  and  earlier  workers)  clearly  separates  the  Cau- 
casus earthquakes  from  the  Kazakh  and  southwest  (JSSR  explosion  population. 

This  line  would  also  separate  the  suite  of  Asian  events  studied  by  the  Dahlman 
et  al.  (1974).  The  Kazakh  and  southwest  USSR  explosions  are  located  approxi- 
mately 10  to  25  degrees  north  and  northeast,  respectively,  of  the  Caucasus 
earthquakes. 

Corner  Frequency  and  Long-Period  Spectral  Level 

Figure  15  is  a plot  of  | Qq [ versus  corner  frequency  for  the  Caucasus 
earthquakes  and  Kazakh  and  southwest  USSR  explosions  recorded  at  LASA  and 
NORSAR.  No  long-period  P phases  were  recorded  at  the  long-period  arrays  due 
to  the  small  magnitudes  of  the  events,  so  the  values  of  | | were  estimated 
from  the  spectral  data  of  the  short-period  vertical  component.  Hanks  and 
Thatcher  (1972)  have  shown  that,  for  a given  long-period  level,  Amchitka  explo- 
sions have  a higher  corner  frequency  than  Aleutian  earthquakes.  The  physical 
basis  of  this  discriminant  is  the  smaller  source  time  dimension  of  the  explo- 
sion for  a given  long-period  level.  However,  Figure  15  shows  no  separation 
between  the  Caucasus  earthquakes  and  the  Kazakh  and  Southwest  USSR  explosions. 
There  is  a slight  trend  toward  higher  corner  frequencies  for  lower  |Gq|  for 
the  explosions  and  the  earthquakes.  The  lack  of  a clear  discriminant  could  be 
due  to  characteristics  of  the  sources  examined  in  this  study.  Another  possible 
explanation  is  that,  due  to  the  low  magnitude  of  the  events  studied  here,  a 
clear  asymptotic  value  of  ) Qq | was  not  possible  in  most  cases. 


Blandford,  R.  R. , and  D.  Clark,  1975,  Variability  of  seismic  waveforms  record- 
ed at  LASA  from  small  subregion  of  Kamchatka,  SDAC-TR-75-12  Teledyne  Geo- 
.tech,  Alexandria,  Virginia. 

Dahlman,  0.,  I.  H.  Israelson,  A.  Austegard,  and  G.  Hornstrom,  1974,  Defini- 
tion and  identification  of  seismic  events  in  the  USSR,  Bull.  Seism.  Soc , ’ 

Am. , 64,  607. 
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Figure  15.  Long-period  spectral  level  vs.  corner  frequency  for  Caucasus 

earthquakes  and  Kazakh  and  southwest  USSR  explosions  from  LASA 
and  NORSAR  P recordings. 


Long-Period  Body-wave  Excitation 


There  were  too  few  long-period  P observations  to  determine  average  event 
ratios  of  these  quantities  to  other  phases.  All  of  the  long-period  P observa- 
tions were  close  to  the  noise  level,  as  would  be  expected  for  these  small-mag- 
nitude events.  No  clear  long-period  S waves  were  observed  on  the  recordings 
analyzed  for  the  Caucasus  earthquakes. 

Depth  of  Focus 

All  of  the  Caucasus  events  studied  here,  except  event  6,  had  apparent  pP 
phases  recorded  at  either  LASA  or  NORSAR  or  both  arrays.  Some  of  the  picks 
made  here  are  tenuous  of  course.  No  LASA  or  NORSAR  short-period  data  was 
available  for  event  6;  however,  pP  for  event  6 was  observed  at  a few  of  the 
WWSSN  stations.  The  pP  phases  recorded  at  the  WWSSN  stations  agreed  with  the 
depths  determined  at  LASA  and  NORSAR  for  all  cases  and  in  fact  helped  to 
establish  the  LASA  or  NORSAR  pP  picks  in  some  cases.  Thus  many  low  M -m 
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events  in  this  region  could  be  identified  as  earthquakes  on  the  basis  of  pP 
observations.  pP  was  of  course  not  identified  for  the  explosions,  but  this 
fact  alone  cannot  function  as  a positive  discriminant.  Although  for  crustal 
depth  events  there  is  little  pP  moveout,  when  both  LASA  and  NORSAR  pP  observa- 
tions were  made  for  the  same  event  they  showed  the  correct  sense  of  moveout. 

Complexities 

We  have  computed  the  "complexity"  parameter  as  seen  at  LASA  and  NORSAR 

for  the  Caucasus  earthquakes  and  Kazakh  and  southwest  USSR  explosions  in  the 

manner  given  by  Lambert  et  al.  (1969).  Figure  16  shows  the  complexity  values 

versus  m,  . We  reiterate  that  LASA  and  NORSAR  beams  show  apparent  pP  signals 
b 

in  most  cases  since  this  presense  will  tend  to  enhance  complexity.  The  com- 
plexity numbers  for  the  Caucasus  events  were  generally  higher  than,  but  over- 
lap with,  the  numbers  for  the  explosions.  The  lowest  complexity  numbers  for 
the  Caucasus  events  for  event  5,  a possible  subcrustal  event  in  the  Caspian 
Sea.  The  low  complexity  value  for  this  event  is  likely  due  to  the  late  arri- 
val of  the  pP  signal,  beyond  the  30-sec  signal  window  used. 

Lambert,  D.,  D.  von  Seggern,  S.  Alexander,  and  G.  Galat,  1969.  The  LONGSHOT 
Experiment,  Volume  II.  Comprehensive  Analysis,  SDL  Report  No.  234,  Tele- 
dyne Geotech,  Alexandria,  Virginia. 


COMPLEXITY 


Figure  16.  Complexity  vs.  m.  for  Caucasus  earthquakes  and  Kazakh  and  South- 
west USSR  explosions  from  1.ASA  and  NORSAR  F recordings. 
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Spectral  Raring 


Spectral  ratios  have  been  calculated  at  LASA  and  NORSAR  accord! 
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Figure  18. 


Short-period  P spectral  ratio  corrected  for  attenuation  vs.  irq,  for 
Caucasus  earthquakes  and  Kazakh  and  southwest  USSR  explosions  re- 
corded at  LASA  and  NORSAR.  The  t*  values  are  for  uj-2  source  models 
for  earthquakes  and  explosions  and  id~3  source  models  for  earthquakes 


S/p  Excitation 

As  would  be  expected  for  these  small  magnitude  events,  there  were  no 
short-period  S observations  from  the  Caucasus  events  to  determine  the  ratio 
S/1  ground  disolacement.  No  clear  long-period  S phases  were  seen  either. 


SUMMARY 


Eight  earthquakes  from  the  Caucasus  and  eight  Kazakh  and  southwest  USSR 
explosions  were  examined  in  a seismic  discrimination  context.  Most  of  the 
Caucasus  events  lie  close  to  major  faults  observed  in  satellite  photographs. 

It  was  impossible  on  the  basis  of  first  motions  to  determine  individual 
fault-plane  solutions  due  to  the  small  magnitudes  of  the  events,  but  thrust 
faulting  may  be  indicated  by  the  composite  pattern.  We  were  also  unable  to 
determine  a radiation  pattern  from  the  LR  amplitudes,  probably  due  to  varying 
propagation  effects.  The  following  characteristics  of  the  Caucasus  events 
indicate  that  all  the  events  are  earthquakes: 

- All  the  presumed  earthquakes  had  apparent  pP  phases  recorded  at  LASA, 
NORSAR,  or  some,  of  the  WWSSN  stations. 

- All  the  Caucasus  events  fall  above  the  line  M = m,  - 1.5.  The  Kazakh 

s b 

and  the  southwest  USSR  explosions  lie  below  the  line. 

Recall  that  the  Caucasus  earthquakes  were  chosen  on  the  basis  that  they 
had  the  lowest  M values,  relative  to  m , in  a preliminary  survey.  We  note 

S D 

that  shot  arrays  could  be  used  to  raise  the  M of  an  explosion  by  about  0.3 
M _ units  or  more,  placing  it  close  to  the  earthquake  population  studied  here. 

Complexity,  moment  vs.  corner  frequency,  and  spectral  ratio  were  not 

particularly  useful  discriminants  in  this  data  set.  The  low  magnitudes  of 

the  earthquakes  studied  here  (mb  from  4.4b  to  5.13  and  M from  3.71  to  4.78) 

made  clear  discrimination  difficult  with  the  available  data,  except  for  M -m 

s b 

data,  and  probably  represent  nearly  the  lowest  threshold  of  multi-discriminant 
analysis  for  the  source  region  and  stations  used  in  this  report.  The  instal- 
lation of  high-quality  SRO  stations  in  Asia  will  permit  the  examination  of 
lower  magnitude  events  with  multistation  short-period  data. 


REFERENCES 


Beliayevsky,  N.,  A.  Borisov,  I.  Valvovsky,  and  Y.  Schukin,  1968,  Transcontin- 
ental crustal  sections  of  the  U.S.S.R,  and  adjacent  areas,  Canadian  Jour- 
nal of  Karth  Sciences,  l‘)67. 

Ben-Menahem,  A.,  S.  Smith,  and  T,  Teng,  1965.  A procedure  for  source  studies 
from  spectrums  of  long-period  seismic  body  waves,  Bu 1 1 . Se ism.  Soe . Ame r . , 
55,  203. 

Blandford,  R.  R.,  and  D.  Clark,  1975,  Variability  of  seismic  waveforms  record- 
ed at  1.ASA  from  small  subregion  of  Kamchatka,  SDAC-TR-75-12,  Teledyne  Geo- 
tech,  Alexandria,  Virginia. 

Borusk,  A.,  and  V.  Shalpo,  1976,  Correlation  of  endogenous  processes  in  the 
Great  Caucasus  and  adjacent  part  of  the  Scythian  Platform,  in  Geodyna- 
mics : Progress  and  Prospects,  ed.  C.  Drake,  American  Geophysical  Union. 

Byus,  E.  I.,  A.  I).  Tskhakaya,  and  M.  M.  Rubinshtein,  1968.  Regional  seismic 
zoning  of  the  USSR:  Georgia,  in  Seismic  Zoning  of  the  USSR,  edited  by 
S.  V.  Medvedev,  English  translation,  Keter  Publishing  House,  Jerusalem, 
Israel . 

Dahl  man,  l'.,  I.  H.  Israelson,  A.  Austegard,  and  G.  llornstrom,  1974,  Defini- 
tion and  identification  of  seismic  events  in  the  USSR,  Bull.  Seism.  Soc , 
Am. , 64 , 607. 

Geller,  R.  J.,  1976.  Scaling  relations  for  earthquake  source  parameters  and 
magnitude,  Bu 1 1 . Seism.  Soc . Am . , 66,  1501. 

Hanks,  T. , and  W.  Thatcher,  1972.  A graphical  representation  of  seismic  source 
parameters,  J^  Geophys , Res. , 77 , 4393. 

Ke i 1 is-Borok , V.  i.,  19t>0.  Investigation  of  the  Mechanism  of  Earthquakes, 
English  translation,  American  Geophysical  Union,  Washington,  D.  C. 

Lacoss,  K.,  1969.  A large-population  LASA  discrimination  experiment.  Technical 
Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 

Lambert,  I).,  D.  von  Seggern,  S.  Alexander,  and  G.  Galat,  19b9.  The  LONG SHOT 

Experiment,  Volume  II.  Comprehensive  Analysis,  SDL  Report  No.  234,  Tele- 
dyne Geotech,  Alexandria,  Virginia. 

Neprochnov,  Y.,  1968,  Structure  of  the  earth's  crust  of  epicontinental  seas: 
Caspian,  Black,  and  Mediterranean,  Canadian  Journal  of  Karth  Sciences, 

5,  1037. 

Nowroozi,  A.,  1971,  Seismo-tectonlcs  of  the  Persian  Plateau,  Eastern  Turkey, 
Caucasus,  and  Hindu-Kush  Regions,  Bu 1 1 . Seism.  Soc . Am. , 61,  317. 


REFERENCES  (Continued) 

R1"8d^jJ:,S^:6*a^^^llh°°d  e8t‘«“  »f  *™»t  magnitude.  Bull, 
seismic  events^ 

for  underground  nuclear^xp^si1”^'  ofophys.  mi  SW!Ctra 


